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We consider fluctuations in a perfect irrotational fluid coupled to gravity in an Einstein static 
universe background. We show that the homogeneous linear perturbations of the scalar and metric 
, fluctuations in the Einstein static universe must be present if the second order constraint equations 

■ are to be integrable. I.e., the 'linearization stability' constraint forces the presence of these homo- 

' geneous modes. Since these linear homogeneous scalar modes are well known to be exponentially 

\ unstable, the tactic of neglecting these modes to create a long-lived, almost Einstein universe does 

^ , not work, even if all higher order (L > 1) modes are dynamically stable. 
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I. INTRODUCTION 



' Recently Barrow et al re-examined the stability of the Einstein static spacetime against arbitrary linear 

'"^ fluctuations in the metric and matter and found, surprisingly, that some modes are stable. More precisely, 

they found that given a sufficiently large speed of sound in the background, all non-gauge inhomogeneous 
, scalar modes were neutrally stable (i.e. the fluctuations are not damped), and furthermore vector and tensor 

OO ' modes were neutrally stable on all scales irrespective of the equation of state in the background. The essence 

\ of this stability arises from the spatial compactness of the Einstein static spacetime, i.e. there exists a 

I maximum physical wavelength in this closed space and furthermore the Jeans length is a significant fraction 

of this maximum scale. It turns out that for specific equations of state in the background matter all physical 
O . wavelengths fall below the Jeans wavelength and the modes are thus stable for the usual reasons. 

This surprising non-Newtonian stability for a large class of fluctuations, which was pointed out earlier in 
different and more restricted contexts in p| and Q , is one of the key elements of support for the Emergent 
Universe models proposed recently by Ellis et al in These models explicitly construct spatially closed, 
positively curved, cosmologies which do not bounce and in which inflation (triggered by precisely the famous 
homogeneous instability of Einstein static) is not preceded by an era of deceleration, by contrast to deSitter 
and most other models of closed inflation. Such constructions have no initial singularity, undergo the usual 
inflationary period ending in the usual reheating era, and immediately solve the horizon problem owing to the 
staticity of the initial state. However another key feature of these models is that a finite number of e-foldings 
of inflation occurs over an infinite time in the past. It is critical for this scenario that the homogeneous 
mode, which is always exponentially unstable in linear order, be excluded with high accuracy. 

In the present paper we investigate the initial value problem for second order inhomogeneous fluctuations 
about the Einstein static solution, as made precise below. For some time now it has been known (see, e.g., |3]) 
that the initial value equations of gravity have peculiar properties in a closed universe with symmetries. The 
crucial idea is that there are certain flux integrals, associated with those symmetries, which are constrained 
to be zero. Thus not all solutions to the linear perturbations equations are actually perturbations of some 
exact solution to the Einstein equations. Technically, the tangent space defined by the linearized initial 
value constraints is actually larger than the manifold of solutions. Only those solutions to the linearized 
equations which obey an additional second order constraint can be the linearization of an actual family of 
solutions to the full Einstein equations. This feature of the perturbations of the Einstein equations is called 
"linearization stabilty" . It does not refer to the dynamic stability of pertubative solutions, but to the above 
"stability" of the linearized solutions as genuine approximations to exact solutions of the Einstein equations. 
The nature of these higher order restrictions on the linear fluctuations about a certain background space, 
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which are necessary conditions to ensure the so-called 'linearization stability' of that background space ( see 
), is such that for every Killing vector in the background there is an associated higher order constraint. 
In this paper we focus on the higher order constraint associated with the timelike Killing field of the Einstein 
static initial state. 

We find, in the case of a general irrotational perfect fluid with cosmological constant, that there are no 
nontrivial solutions to both the linearized constraints and this global nonlinear constraint when we exclude 
the linear homogeneous scalar metric and matter fluctuations. In other words, the leading order linear 
metric and matter fluctuations must be trivial if their linear seeds do not include 'zero mode' homogeneous 
fluctuations. It is well known that these homogeneous modes are dynamically unstable for any (causal) 
equation state of the matter, just as all perturbations of an Einstein static dust model are unstable. This 
would seem to suggest that if the universe is in a neighbourhood of the Einstein static solution then it does 
not stay there, even if perturbed with the neutrally stable modes found in reference H. This is one of the 
first instances where the linearization stability issues make a physical difference in the analysis of metric and 
matter stability (see also Brill, 

The paper is organized as follows. In section II we briefly outline the details of the Einstein static 
background model. In section III we define the constraint equations in the standard ADM decomposition 
of the Einstein equations and also define an orthogonal decomposition (following Q) of perturbations into 
transverse and longitudinal parts. In Section IV we formulate and compute the nonlinear constraints, leaving 
conclusions for Section V. The entire analysis is quite similar to that in Brill and Deser's original paper Q . 



II. EINSTEIN STATIC SPACETIME 



Consider a FRW universe in comoving coordinates {t,x) with scale factor a{t), with signature (-1,1,1,1), 
and with a perfect fluid with energy density p and pressure p. The equations of motion for the scale factor 
a(t) are, according to the Einstein equations, 

- = -|b(l + 3«;)-A], (1) 
a 6 
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= -ip + A)^-, (2) 



where K = ±1,0 is the constant curvature of the t = const slices, H = dtln{a) is the Hubble parameter, A 
is a cosmological constant, w = ^, and k = SnG in units where c — 1. One can combine these equations to 

fovmH = -q{p{l + w)) + §. 

As Einstein found some time ago, demanding that the universe be static (a = a = ) obviously sets K to 
be positive (take it to be 1) and leads to constraints relating the initial energy density and pressure of the 
fluid to A. The equilibrium radius of such a static universe is set by these constraints to be 

' (3) 



" A(l + w)' 

where (1 -I- w)p > 0. In the case of dust it is seen that Oq = 1/A. 

Consider the general perfect fluid equations. The equations of motion of the perfect fluid are given by the 
conservation equations of T'^'' — [p + p)u'^u'^ + pg'^'^, or 

{p + p),^,u^'u^ + {p + pX^"!. + (/O + P)u''u^;i, + P;,. = (4) 

Dotting this with u'^ we get the primary conservation law 

Multiplying this by Ui, and subtracting from the above equation we get 

= {p + p)u'^u^-f_, + p-i, + p^^u'^Ui, 

= {P^)iVNu[^).,,^u'^, (6) 
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where 

N{p) = exp (^2 1 (7) 

is a function of p where we assume that there exists some constitutive equation p — p{p) . These equations 
are trivially satisfied if we assume that the term N^{p)u^ satisfies 

n"^{p)u^^(^,^,, (8) 

where is a scalar potential for the flow. This is an appropriate generalisation of the irrotational conditions 
of non-relativistic flow. Finally, multiplying this equation by itself we have that 

N{p) = (9) 

which is the equivalent of Bernouli's equation for non-relativistic irrotational flow. We will concern ourselves 
with constitutive equations where p is linearly proportional to p. 

We wish to treat the case of an arbitrary irrotational perfect fluid in terms of this scalar field <j) that acts 
as the velocity potential of the fluid (as in 9] ). We choose the scalar field to have the action 

5 = ^ y(-0'V,Ary^rfS, ae5R (10) 

Comparing the stress energy that results from varying the above action with respect to the metric with that 
of a perfect fluid, we can easily identify the corresponding energy density and pressure in terms of 0: 

P = io^-\)N^. (11) 

V - - (12) 

For a = 1 we obtain the stiff {w — 1) perfect fluid, i.e. a minimally coupled scalar fleld. However, note that 
the speed of sound is given by = p/ p = (2a — 1)^^ and thus causality restricts a > 1. It turns out that 
stable (inhomogeneous) scalar fluctuation modes only exist when the background speed of sound satisfies the 
bound v1 > ^ (^, 0) for no Jeans instability, which translates into a < 3. I.e. there is no Jeans instability 
for any physical inhomogeneous modes given these conditions. Thus we take 1 < a < 3 in what follows. We 
also note in passing that in order for the four velocity to be timelike the gradients of the velocity potential (j) 
are restricted to ones such that their temporal gradient dominates their spatial counterpart. In this model 
we take the background — ( ^^°^^ )~t to satisfy equation (7) in the background, so the timelike condition 
will always be satisfied for the perturbations. 

In any case the generalized Einstein static initial conditions (in terms of (f) become 

^ a + 1 . a + 1 . ^ 

TV"—— = A = ^, 13 

2 anal ^ ' 

i.e. = (a + l)/(aA) and N = (2A/(K(a + 



III. INITIAL- VALUE CONSTRAINT EQUATIONS 



The gravitational field may be characterized in terms of the initial three-geometry and extrinsic cur- 
vature K^^ of some spacelike surface. These twelve quantities may be rewritten into a more convenient 
hamiltonian form by defining the momentum density tt*-' = ^J\g\(K'^^ — K^^g^^), which is conjugate to gtj. 
The phase space for Einstein's equations is then some suitable function space of pairs [gab, t^""^) over a three- 
dimensional manifold M (which we take as compact and without boundary), and the constraint subset of 
phase space is defined by the four initial-value constraints. The constraints are nonlinear and we expand 
them order by order in the strength of the metric and matter perturbations. Quantities with overbars will 
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indicate zeroth order, background, quantities and a bar indicates covariant differentiation with respect to 
the bacfcground metric (used intechangably with the standard Vi notation). 
The constraints with matter sources (p, J'^) are 

^ ^TT^'nab - + = 2«;/Mp, (14) 
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in-") |fc = «J^ (15) 



where R is the Ricci curvature scalar for the three surface (in what follows we will drop the 3 superscript) 
and TT = gabT^""^- 

At zeroth order, i.e. in the background Einstein static spacetime, there is only one nontrivial constraint 
namely 

3A , , 

where we have used the background conditions Rab = A^ab, tt"'' = 0. At linear order (using the background 
equations) the constraint equations are 

5R = 2k6p, (17) 
VbtSTT'^'' = n^\5r, (18) 

and similarly the second order equations are, using the linear and background equations, 

= S-K-'^STTab - IsttA + 2k(5V (19) 

m V 2 y 

VbS^TT'^' + 2 [SCU^^^' + '5CL<57r-] y/Wl = ^S'ir^/\9\), (20) 

Here ^SC^^^ = g°''^{y mSgbc + '^bSg-mc — c5gbm) is the perturbed connection. One can rewrite equation (19), 
the equation associated with time reparametrization invariance at second order, using 5'^R = 26g'^''SRab + 
g-'^S'^Rab + S^g-'^Rab to obtain 

K5^g - ^"'y'S^gbm - AgabS^g-' = ^ iB' + 25Cl^JCl^,r'' + ^dg'^'SRab - (^Sn'^'STTab - ^Sn^^ - 2nS^p, 

(21) 

where by 'SI ^B^ we denote a combination of terms that occur in the form of a total derivative. In obtaining 
(21) we have expanded the second order Ricci scalar in terms of the perturbed connection defined above and 
grouped second order terms and products of linear order terms. One can similarly simplify the momentum 
constraints (20) but in this paper our main concern will be the second order timelike, or "Hamiltonian" , 
constraint. 

Typically what one does is to solve the constraints for a given set of fluctuation variables and in that way 
solve for the constrained variables in terms of the free variables. A convenient way to facilitate this process 
is through an orthogonal decomposition of the fluctuations in gab and n""^ into longitudinal and transverse 
parts. As pointed out in Isl] this procedure differs crucially for tensors in closed spaces from the usual one 
in flat space since one can now have fluctuations in global parameters such as the total volume of the space 
while still having no local excitations. 

Following the notation of Brill et al one may decompose the metric perturbation into its various transverse, 
transverse traceless, and longitudinal parts via 

5gab = 5gZ^^ + \ (-VaVb + K-gab) -Jg^^''' + A&fgab + 25g[^)), (22) 

where A = VcV"^. The transvserse traceless and longitudinal parts are defined as usual: 

= = V''55ir\ (23) 



V,<55(^)" = 0, (24) 
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and here Asg = Asg + ^dg^^^\ The homogeneous fluctuations Asg are essentially the global fluctuations 
(related to the volume fluctuation) unique to perturbations of closed spaces. An analagous decomposition 
holds for the momentum fluctuation (the homogeneous modes Asn describe the 'time rate of change' of the 
volume fluctuation). 

Using the above decomposition we note that an identity due to the symmetry of the background is that 
the perturbations are all 'doubly transverse'. Indeed, we see 



(A + A)(5<?i' 



We also note that we can write 



{V) 



S'R = AgabS^g'^' + 25g''''5Rab + ir'^ClyJC\, - 2^^J'Cyg'^'> + V,B' 



where 



25^C\, = {Vb5^gcm + yc5^gbm - ym^^gbc) 

We may use the decompositions above to eliminate the 'doubly transverse' term and obtain 



(25) 
(26) 

(27) 

(28) 



1 



1 



(29) 



where 5^ p are the second order energy density fluctuations. We explicitly insert the above decompositions 
into equation (28) in the next section. 

In the case of the generalized scalar fleld (f) introduced in section II, the Hamiltonian formulation is more 
involved since the canonical momentum conjugate to is generally 



dC 



d(j) AA 



^(.^-A^V.i) 



A/'2 



2A/'^ 



ct-l 



A/'2 



where Af and A/"* arc the purely kincmatical lapse and shift variables of the ADM formalism. For general 
3 > a > 1 this expression cannot be inverted in closed form to find <p = (t>{'J^4>), however for a spatially 
homogeneous scalar fleld it is possible and the resulting Hamiltonian (energy density) is 



Ho 



2a - 1 / {tt^Y 



2a 



A 
+ - 



(30) 



l2q-1 



where tt^ = -v/l^li?!! . We define the fiuctuations Jtt^ and (J^tt^ using the formal definition of -k^ given 
above. For example the first order fiuctuation in is formally defined by 



2a- 1 



2{a - l)5N 



(31) 



and similarly for the second order fluctuation (though new cross-terms like SAf'^Scfy^i and S^'''5cj)^i start to 
appear). Formally perturbing equations (11) and (12) with the above deflnitions of dir^, S'^tt^ yields 



2k6p 
2kS'^P = aA 



, 6N 2Aa 



(a-1) 



SN 



Sw^ _ Sg^ 

7f<Y, 2 

+ 



N 



(32) 
(33) 
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where the second order term is given expUcitly by 
so that, in Hamiltonian form, we finally obtain 



2aA 
2a -1 



\/W\ 



+ j2^^ZrY)^ ((2a2 + 8a - 7)((51n ^/Wif + {2a^ + 4a - 5)((51n7r0)2 - A{a^ + 3a - 3)((51n ^f\g\){&\n 

Here 5g = g°'^5gab and TV" = A/(K(2a — 1)) by the zeroth order constraint (13). Inserting the above matter 
perturbations into equation (25) and using the second-order identity g°'''S^gab + gabS^g"''' — —26g°-^5gab^ we 
arrive at 



(A + A^r^)5'g = A il-^ ) Sg'^'Sgab + 25g'^''5Rab + 2<5C,%<5C,%g 



2a- 1- 



2a- 1 



2Aq! 
'2a- 1 



7f0 



1 



(2a -1) 



31 V 2 
(-2a2 + 10a-8)((51nA/|ff|)^ 

(35) 



+ {2a^ + 4a - 5){SlTnr^f - ^{c? + 3a - 3)(^ In y/\^\)(i\^T:^)^ 

The various terms in ((Jtt^) can be reexpressed in terms of the metric fluctuations, using the linearized 
Hamiltonian constraint, via 



(5 In TT^ = 



(2q-1) 
2Aa 



A + A 



a — 1 



2a- 1 



bg - ^'^'Sgcb = MSg 



(36) 



so that the only explicit matter dependence in the second order hamiltonian constraint appears through the 
S'^TT^ and g^^ 5(j)^i54>^j terms. 



IV. NONLINEAR RESTRICTIONS ON THE LINEAR FLUCTUATIONS 

A. Fixing the linear gauge 

As is usual in relativistic perturbation theory one must fix the coordinate freedom inherent in the metric 
and matter fluctuations. In our particular case we want to remove the homogeneous second order matter 
dependence, which enters through the term, in order that equation (35) is of the form C5^F = S{{SF)'^) 
(where C is an elliptic operator (with only constants in its kernel) acting on the fluctuations in some quantity 

in a closed space). In this form (since the timelike Killing vector component is trivial) the nonlinear 
constraints simply become J S = where the integral is over the closed space @ . 

Fortunately the Einstein static background crucially simplifies the relevant gauge transformation laws of 
the fiuctuations not only at linear order, but also at second order. Indeed the term (5^7r^ only depends on 
the linear gauge fixing essentially because tt^ = 0. We are only interested in the homogeneous part of this 
S'^TT^ term so we pick a linear gauge-fixing to eliminate it and leave only the inhomogeneous part S'^tt^. 

We consider the linear spacetime coordinate transformation — + (T, + M*), where (T, M, M') 
satisfy 



A((2a-l)A + (a-l)A)r = -{2a - 1)Stt^^'''> 
I {2MASTr^+S^TrA ^/\f\(fx = 



(37) 
(38) 
(39) 
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and where the spatial dependence of the modes is understood in terms of the eigenfunctions of the spatial 
laplacian: AF ^ -J^l ^^^"^^^ = ^^(^ + 2)-^, Le Z+. In this new coordinate system Sg'^y^ = 0, 

the homogeneous part of S'^tt^ is zero, and 5t:'^^^'^ — 0. Thus 

^ (k afcc 1^ 2\ 2Aa 2 1 fx abc TT 3 ~ 3^ ^Ag 2- /.n^ 

- — Uvr JvTab - -<57r - — —6 tt^ = dnTj^dn^^ - - Att) - ^-7: tt^. 40 

\g\\ 27 7r0(2a--l) |g| V 2 J Tr^{2a ~ 1) 

This remarkable tranformation completely fixes the vector degrees of freedom at linear order, however we still 
have the freedom M M + f(t) in the scalar sector which we use to eliminate the homogeneous modes Asg 
by picking a special f{t) (whose form is not particular illuminating at this stage). Furthermore, the linearized 
Hamiltonian constraint (36) in this coordinate system implies that AJtt^ cannot be zero everywhere, which 
means we can always pick a function M such that equation (38) is satisfied. 



B. Gauge-fixed nonlinear constraints 

The nonlinear constraint associated with (35) (an integrability condition on is effectively the 

integral of the right hand side of equation (35) set to zero: 



= 



A 



S3 



2 -3a 

2Q!- 1 

2Aa 



2a -I 



+ {2a' +Aa- 5)(J Iutt^)^ - 4(a^ + 3a - 3)(51n Vl5|)(<51n7r0)JJ y/\g\d 
Using by-parts integration and compactness one can show that 



-g\[5g'''6Rab + SCt^JCi^,r') = j ^Jg^'SRab^, 



which implies, using the linear equations of motion and making the above gauge choice. 



1 



-TT , x„ab^^9ab 



Snl^PrSn'J + Sg" 
S3 IffI 2 



A 



2(2a - 1 



■Sg'^'Sgab 



2Kc 



(2a- 1)3 



M^{2a'^ + 4a - 5) - 2M{a^ + 3a - 3) + 



-2a2 + 10a 



{5gf 



(41) 



where V represents the volume of the initial static space and Ai is defined in equation (31). We have 
eliminated the terms in 5(l)^i by using the scalar-sector momentum constraints (which are ^ ( Stt'^'^'^^ — 

S(j)) = ). The constraint is split into its inhomogcncous {L > 2) and homogeneous {L = 0) pieces. [To| . 

Finally, inserting the decompositions into equation (41) and using equation (36) to remove the dependence 
in (Jtt^)^, we arrive at 



= / A^/\f\d^x-^3 



A' 



2 



V, 



(42) 



where 

E 

L>2 
+ 



3 2 / a^(2a2 +4a - 5) 6 a(6a2 + 21a - 20) 4 lOOa^ - 194a2 + 109a + Sa^ - 20 



aA'^fc 
a- 1 



2(a-f l)4(2a- 1) 4(2a - l)(a + 1)3 

1 



(a+ l)3(2a- 1)3 



(43) 



2(2a- l)(a + 1) 



(tr)N2 



\9\ 



+ -^ia{e -2) + ^±l)Sg<,^^^Sg^r> ^ g^'^g 



ia -mb 
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and where = L{L + 2). We point out that equation (42) is not the integral of the second order Hamiltonian 
for the fluctuations since, given our gauge fixing, the symplectic terms Sir^'^^^dgl^^ do not contribute. 
The main result of this paper is that first term in equation (43), A, is positive definite given ct > 1, L > 2. 
We observe that the global constraint is an integral over S^, i.e. is of positive measure. This means that in 
the absence of the homogeneous mode, which provide a strictly negative definite contribution to the integral 
through the second term, there is no nontrivial solution to these global constraints even though there are 
certainly solutions to the second order equations which only have inhomogeneous linear seeds. In other 
words, if one wants to study the evolution of the second order modes one must include, as part of their 
source, the zero mode at linear order in order to properly satisfy the initial value constraints. Therefore, we 
have shown that at linear order we must include the unstable, homogeneous mode. We emphasize that we 
are not claiming stability or not at the second order level, rather we are claiming instability at the linear 
level. Whether or not second order perturbations could stabilize the spacetime at a sufficiently high value 
of the linearized solution remains unclear, but we do not address that issue here though we strongly suspect 
not. 



V. CONCLUSIONS 

We have shown that the requirement that the second order Einstein constraint equations be integrable 
demands that any non-homogenous linear mode perturbations of the einstein static universe must be accom- 
panied by the homogenous linear mode with comparable amplitude. Since this homogeneous linear mode 
is exponentially unstable, this implies that any linear approximation to a solution of Einstein's equations 
must be unstable. Our result is valid for perfect fluid matter determined by a potential and with a constant 
velocity of sound. 
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